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Abstract We investigate the diffusion of particles in an attractive one-dimensional 
potential that grows logarithmically for large \x\ using the Fokker-Planck equation. 
An eigenfunction expansion shows that the Boltzmann equilibrium density does not 
fully describe the long time limit of this problem. Instead this limit is character- 
ized by an infinite covariant density. This non-normalizable density yields the mean 
square displacement of the particles, which for a certain range of parameters exhibits 
anomalous diffusion. In a symmetric potential with an asymmetric initial condition, 

, K the average position decays anomalously slowly. This problem also has applications 

outside the thermal context, as in the diffusion of the momenta of atoms in optical 
molasses. 
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1 Introduction 

Diffusion of particles in the presence of an attractive logarithmically growing po- 
tential has attracted much interest since it serves as a model for diverse physical 
systems. Examples include charged particles in the vicinity of a long and uniformly 
charged polymer[22| (the Manning condensation problem), diffusive spreading of 
momenta of two-level atoms in optical latticesl 8 , 23 , 2Tl[T2ll . single particle models 
of long ranged interacting systems[6 ,5 ,9], probe particles in a 1-d driven fluid[20'|, 
dynamics of bubbles in double stranded DNA I13lll4l l3l l30l . vortex dynamics Q, 
and nano-particle dynamics in appropriately constructed force fields [ 10 1. When the 
motion is in contact with a heat bath, the steady state is given by the equilibrium 
Boltzmann distribution, W eq (x) « exp[—V(x)/kBT]. For a logarithmically growing 
potential, V(x) ~ Voln(x/a) forx ^> a, Vq > 0, we get an equilibrium distribution with 
a power-law tail, W eq (x) ~ x~ v °' kBT . By definition this steady state is normalizable 
when Z = f™„dx exp[~V (x) / ksT] is finite, a necessary condition being Vo/ksT > 1. 
The well investigated Bessel process [24 27] considers a purely logarithmic poten- 
tial, V(x) — Voln|x|, which is not regular at the origin. This means that the steady 
state is not integrable, diverging either at large or small x, depending on whether 
Vo/ksT is smaller or larger than unity. In most cases, the potential is only asymp- 
totically logarithmic, with a regular small x behavior characterized by some length 
scale a. For example, this length is the diameter of the long polymer in Manning's 
problem. Similarly for optical lattices, the potential (here, in momentum space) is 
V(p) = !/2ln(l +p 2 /pn), so that po sets the length scale. 

Here we study the case of diffusion in such a regular logarithmic type poten- 
tial. We mainly discuss the case where there exists a time independent steady state, 
i.e., the Boltzmann equilibrium distribution is normalizable. Nevertheless, higher mo- 
ments of the distribution are infinite, due to the slow spatial decay of W eq . For exam- 
ple, the second moment, {x 2 )eq = J- aa x 2 W eq (x)dx, diverges for Vo/kgT < 3, and the 
fourth moment for Va/kgT < 5. The equilibrium distribution is not sufficient then to 
compute all moments higher than some given moment, depending on the size of Vq. 
Rather, as we will see, the higher moments grow in time, approaching the equilibrium 
value of infinity in the infinite time limit. At any finite time, however, the moments 
are perfectly finite, as the physics demands. The computation of these moments in 
the long-time limit involves a new object, which we call the infinite covariant den- 
sity. The term infinite in the name |29, 1, 19] refers to the fact that this density is not 
normalizable due to its divergence at the origin. Nevertheless, the higher moments of 
this distribution are well defined. The term covariant refers to the fact that the den- 
sity has a scaling form, which gives rises to the power-law growth of the moments 
with time. The scaling properties of the distribution and the power-law growth of the 
moments have been obtained in Refs. B20II171 by assuming a scaling ansatz. Here we 
construct the time-dependent probability density from first principles via an eigen- 
value expansion of the Fokker-Planck operator. This new derivation will allow us to 
calculate additional properties of the system, in particular the anomalous decay of the 
first moment from an asymmetric initial condition. This long-term memory is a sign 
of the breakdown of ergodicity for too small Vq ifTTI . 
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The plan of the paper is as follows: We first construct the Schrodinger opera- 
tor associated with the Fokker-Planck equation, and express the probability density 
function (PDF) in terms of the eigenstates of this operator. We then specialize to 
the exactly solvable case of a pure logarithmic potential connected to an inner re- 
gion of constant potential. This allows us to explicitly solve for the eigenstates and 
construct the PDF, and to exhibit its scaling properties in the long-time limit. In the 
following section, we generalize the procedure to an arbitrary logarithmic -type po- 
tential, which reproduces our previous scaling solution [ 17 ] for the PDF. We calculate 
all even moments to leading-order, and the first-order corrections for the zeroth and 
second moments. In the final section, we discuss the decay of the first moment. We 
conclude with some final observations. 



2 General Formalism 

We consider Brownian particles in a potential as modeled by the Fokker-Planck equa- 
tion (26) 

jW(x,t)=L w W(x,t). (1) 

Here W(x, t) is the probability density function (PDF) of finding a particle at x at time 

t and 

d 2 d , , 
-F(x) 



L FP =D 



dx 2 dj 



(2) 



subject to some initial density Wq(x) at time t = 0, where D is the diffusion con- 
stant. If the system is close to thermal equilibrium, the Einstein relation gives F(x) = 
f{x)/kBT where T is the temperature, and f(x) is the external force acting on the 
particles. For other processes not necessarily close to thermal equilibrium, e.g. atoms 
in optical lattices 023], the temperature is not defined in its usual sense, and in what 
follows we do not limit ourselves to thermal systems. The force F(x) is derived from 
a potential U(x), F(x) = —U'(x), which for a thermal system is related to the poten- 
tial energy V(x) by U(x) = V(x)/ksT, Then, the normalized steady state solution of 
the Fokker-Planck equation is given as usual by 



i /.» 

Weq (x) = -e- u ^- Z^j 



dxe- u W , (3) 



which for a thermal system is the Boltzmann distribution. 

The formal solution of the Fokker-Planck equation is given in terms of the eigen- 
function expansion |26l 



W(x,t) = e- u ^' 2 £ a kWk {x)e-^' , (4) 

/t=0 

where 

a k = J" dxW Q {x)Wk(x)e u ^ 2 . (5) 
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In particular, for a 5-function initial condition, Wq(x) = 8(x — xq), a^ = e u ^ x °'' 2 \i/i ( (xo). 
We will mostly concern ourselves here with the case xq ~ 0, examining the more gen- 
eral situation at the end, in which case 

W(x,t) = e- u ^ 2+u ^ 2 £ Yft(0)Y*(*)«rV (6) 

k=0 

The wave functions %(.*:) satisfy the Schrodinger equation 

&y k (x) = -X k y k (x) , (7) 

where the effective Hamiltonian is the similarity transform of the Fokker-Planck op- 
erator: 

" d 2 



H = e v ^l 2 l w e- V ^I 2 =D 
so that 



jh2 -i/2F>(x)-i/4F 2 (x) 



(8) 



V eff (x) = D (i/2F'(*) + i/4F 2 (x)) (9) 

is the effective potential of the Schrodinger problem.. 

One eigenstate is immediate, namely the Xk = bound state, which is simply 

VAoW=Z- 1 /V t/ W/ 2 (10) 

so that «o = Z -1 ' 2 , independent of the initial state Wq(x). This corresponds to the 
ground state of the effective Hamiltonian (since it has no zeros) and describes the 
steady state. Indeed the k = term in Eq. (J6l) is the steady state Eq. p). In many sit- 
uations (e.g., a binding harmonic potential) the eigenspectrum of the Fokker-Planck 
equation is discrete and then in the long time limit the steady state is the solution 
of the problem, with corrections that decay exponentially in time. In other cases, the 
spectrum is continuous, but then a steady state is not normalizable (e.g. free particle 
on the infinite line). 

We are interested here in the case of a (reflection symmetric) logarithmic poten- 
tial, defined by its asymptotic large-x behavior, U(x) « Uoln(x/a) +o(l). This case 
is different from the more typical examples discussed above. For Uq > 1 there is a 
continuous spectrum starting at zero, since F, F' go to zero for large x and so the 
effective potential of the Schrodinger equation, given by Eq. d5), vanishes for large 
x. In addition, there is the zero energy ground state, since a normalizable steady state 
exists, lying at the edge of the continuum. As we will demonstrate, the steady state, 
and hence the Boltzmann equilibrium distribution for a thermal system, does not fully 
describe the density profile of the particles at long but finite times in the case of the 
logarithmic potential with Uq > 1 . For Uq < 1 , on the other hand, there is no normal- 
izable steady state, and the spectrum is pure continuous extending from E = + .. For 
a general logarithmic potential, of course, the continuum states cannot be solved for 
exactly, so we proceed in the next section by solving a specific simple example. 
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3 A Solvable Example 

A simple potential which admits an exact solution is 

^Ll i; h; a , (id 

v ' 1 Uom\x\/a \x\ > a 

which for the thermal case is associated with the potential energy V(x) = 9(\x\ — 
a)kBTUoln(\x\/a). The associated F(x) in the Fokker-Planck equation is 

n«> -{_% »<: , 



and the steady-state PDF is 




with 



2aUo 



(13) 



(14) 



Uo-l 

The logarithmic potential yields a density profile which is a power law. As noted 
above, for the steady state to exist we must have Uq > 1, otherwise the steady state is 
not normalizable. Notice that in the steady state the second moment (x 2 ) = J ^ dxx 2 W eq 
°° diverges if 1 < Uq < 3. Physically for any finite long time, with the particles starting 
on the origin, the mean square displacement cannot be larger than what is permitted 
by diffusion (x 2 ) < 2Dt. In what follows, we consider Uq > 1 and initially at time 
f = the particles are located at the origin x = 0. 

For our solvable example, Eq. ( [TT] ), the Schrodinger equation |7]i reads 



D 



d 2 «/, , ■,f u o+ 2u o\ 


5(1*1 —a) , r 


^ +0(W " a) ( Ax 2 ) 


~-^i -U Q 

la 



Wk = hWk, (15) 



where the Heaviside function 9{\x\ — a) = 1 if |*| > a and zero otherwise. As men- 
tioned above, we choose for now an initial condition with the particle at the origin, 
from which it follows, given the symmetry of the potential, that the solution W(x,t) 
also has even parity. Hence in our case it is sufficient to investigate only the even 
solutions of the Schrodinger equation. As we discussed in Section II, the spectrum 
consists of the ground state with zero energy and a continuum that goes down to 
A* = + . We thus write the energy of the continuum states as A^ = Dk 2 , k > 0. 

In the region |x| < a the solution is trivial, since the effective potential in the 
Schrodinger equation is zero, we have the behavior of a free particle, whose even 
solution reads 

Yk(x) = Nk cos kx for |jr| < a. (16) 
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Nk is a normalization constant to be determined later. For x > a, the general solution 
with Xt = Dk 2 is 



with 



\j/ k (x) = N k 4 / - A^y a (foe) + £*./_ a (foe) 



l + I/o 



x > a 



(17) 



(18) 



and since the wave functions are even it is sufficient to consider x > 0. 

The coefficients A* and B^ are determined from two boundary conditions at |x| = 
a. First, the continuity of the wave functions, Eqs. ( [To} and (17) , gives 



A/ l J a (ka) +Bi i J- a (ka) = cos(fca). 



(19) 



The 5 function in the effective Hamiltonian gives a discontinuity in y/' k at x = a, so 
that 



dWk 
dx 



dWk 
dx 



Uo 



Wk(l) 
la 







(20) 



Eqs. ( |19|20| i determine the two unknowns A* and B^, which in principle gives a formal 
but rather cumbersome solution to the problem. Luckily we are interested only in the 
long time limit of the problem where naively one expects the steady state to describe 
the system. For that limit, it is sufficient to consider only the small k — > behavior, 
since for long times the modes with finite k decay like exp(— X^t ) = exp(— Dk 2 t ) and 
hence are negligible. More formally, the summation over the eigenfunctions in Eq. 
(6| must be replaced with an integration since the spectrum is continuous, (see Eq. 
( 28] > below) and hence we may consider Eq. (J6l) as a Laplace transform (where t is the 
Laplace variable) because of the exp(— Dk 2 t) term. Then according to the Tauberian 
theorems, the long time behavior corresponds to the k — >• limit. 



In the small k limit we use Eqs. ( 19|20i and the small z limit of J a {z) E), 



Ja(z) 



(!)' 



r(l + a) 



to find 
where 



Ak ss g/(ka 



\2-a. 



4(1 + a) 



Bi.: 



.^ = -2 a - z r(a-l){2a-l); 



m = T a T{\-a). 



(21) 

(22) 
(23) 



In particular, A k ^> B k for a > 1 . 

Now we determine the normalization N^ from the condition 



1 



[Wk(x)] 2 dx 



{N k 



dx cos kx+ / dx(x/a)[Ai c J a (kx)+B] c J- a (kx)Y 



(24) 



where L — >• °° is the system size (a parameter which will eventually cancel out). The 
second integral grows with L, and so only the large x behavior of the integrand is 
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relevant. Using the large z asymptotics, J a (z) ~ ^/2/(7Tz)cos(z — an/2 — 7T/4), and 
the fact that Ak ^ Bk, we find in the large L limit 

^2^L^ 2a " (25) 



Having found the coefficients Aj, B^ and N% we insert the eigenfunctions Eq. ( 17 1 
in the expansion Eq. (|6|. We separate the time independent part (the ground state 
k = 0) from the contributions of the continuum states (k > 0) and write 

W{x,t)=W eq (x)+W*{x,t) (26) 

where, for x > a, 

W*(XJ) « ( X V /2 - a y * (ka) 2a-3 e -D* x 



W- [^Ofca) 2 - a / a (fcx) + ^(/fca) a /_ a Ofcx)] . (27) 

In the large L limit, the sum of k becomes an integral (requiring % to vanish at x = L, 
for example, implies an asymptotic spacing of Ak — n/L) so that 

fc>0 Jo % 



We can now rescale k in the integral by \/Dt and investigate the time dependence of 
W* . In particular, for |x| < a, one finds that W* decays as t a , and so is smaller than 
W eq in this domain. For \x\ > a, the J a term scales as t~ a i 2 whereas the J- a term 
scales as f 1-3 "' 2 , which decays faster for a > 1, i.e. Uo > 1. Thus, forx > a, using 
Eq. 6.631.5 of Ref.HSl 



w*oc,o « ^x i - a fl 2 "- 2 (Df)- a/2 r dkk a - x j a ( kx 



2srf K ' Jo \WDt, 

= - r ( ,-i)U-i). i^r"^\rw-n«A\ ■&> 

Here 

r 

(30) 



V4Z)f 

is the scaling variable, which is immediately familiar as Brownian scaling, and we 
are considering the limit of large times with z fixed. r(a,x) is the incomplete Gamma 
function|2|. Plugging into Eq. (26 1 and rewriting W eq in terms of a gives 

/ s Ot-1 /X\l-2a -1 /4Df\ 1/2_a , ,„ r . , , ,.-, 

WM * W^TTa (a) + r(a-l)(2«-l) fl UJ Z ^ ^ ^ )] 

* v ' "> v ' 

W«,,W W*(x.r) 

(31) 
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The first term on the right hand side is the ground state, namely the usual Boltzmann 
equilibrium. Crucially, this is exactly cancelled by the first term of W*, which is 
essential since W is cut off and thus much smaller than W eq in the large x region. We 
are then left with 

1 /ADt\ x l 2 - a 



r(a-l)(2a-l)a \ a 2 J 

which clearly differs from the steady state Eq. ( p"3} . For z small, the F function is near 
unity, and W reproduces the (x > a) equilibrium state. For z large, however, r(a,z 2 ) 
decays as a Gaussian, 

r(a,z 2 )«z 2(0 ~ 1) *~'\ (33) 

ensuring the finiteness of all the moments. For large z, then, 

The scaling structure of W, namely W(x,t) = t l ' 2 ~ a ^(z) is exactly the starting point 
of our scaling ansatz presented in Ref. 1 17 1 

We have assumed in the above that a > 1, i.e., Uq > 1. The situation for a < 1 
is markedly different. Here there is no normalizable equilibrium solution. Hence, the 
k = mode is not a part of the spectrum of the Fokker-Planck operator. Also, B^ ^$> A^ 
in the small-£ limit, and so 

{Nk)2 ~ 2WL {ka)l ~ 2a ■ (35) 

Thus, for a < 1, we have, forx > a, using Eq. 6.631.4 of Ref. 031 . 



1 ,_,„/4DA 



r{\-a)a 



-1/2 



l - 2a I ™ e~* . (36) 



For fixed x, so that z<^ 1, it behaves as t~ l+a x l ~ 2a . In other words, the central region 
behaves like the nonnormalizable equilibrium solution, with a prefactor that decays 
in time. We also see this by looking at the behavior for x < a, where 

1 «*>r , ,37) 



r(l — a)a \ a 2 J 

which, since W eq (x) is constant for x < a, is again the equilibrium solution with a 
prefactor that decays in time. 
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4 General Logarithmic Potential 

We now turn to the general case of a logarithmically growing potential, which will be 
seen to have the same basic structure as the above calculation. Consider a symmetric 
potential U(x) in such a way that as \x\ — > °° we have U(x) ~ t/oln(|x|/a) + o(l). 
We set (7(0) = 0, since a constant shift in the potential can be accommodated by a 
shift in a. We also assume a steady state exists, that is the potential is sufficiently 
regular at the origin that Z is finite and so W eq (x) ~ Z (x/a)~ u ° for x — > °o, The 
partition function Z depends on the exact shape of the potential, and it is found from 
the normalization condition of the steady state: Z = 2 J °° e~ u ^ x 'dx. 

Again the key to the calculation is finding the continuum eigenstates for small 
k. For large x, we can approximate the potential by its asymptotic form U(x) sa 
U()ln(x/a), and so the Schrodinger equation is the same as we solved above, giving 

Wk{x) =N k J- [A k J a {kx) +B k J- a {kx)} , (38) 

V a 

where again a — (I +Uq)/2. This approximation clearly fails for the central region 
where x is of order a. Instead, we will solve the Schrodinger Eq. as a power-series 
in k 2 . This series breaks down, as we will see, at large x. Nevertheless there is an 
overlap region n«i« \/k where both approximations are valid and will enable 
us to determine A k and B k . To leading order we can drop the Dk 2 \y k term in the 
Schrodinger equation, so that Hl//j w 0, and we have that y/ k is none other than the 
zero-energy solution y/b, up to normalization, 

Wk ^N k e- u ^' 2 . (39) 



For large x, this behaves as \jf k w N k (x/a) <- a , which comparing to Eq. (38 i for 



small kx gives us, using the small argument expansion of J± a , Eq. (21 1 



B k ^m{kdf\ 38 = 2- a r{l-a) . (40) 

This agrees with what we found in our special case in the previous section. To deter- 
mine A k , however, we need to go to next order in perturbation theory in k 2 . 
Writing 

Wk(x) = Wk,oi x ) + (ka) 2 Wk,i{x) + ■ ■ ■ , (41) 

where \j/ k o (x) = N k e~ u ( x >' 2 is the zeroth order solution. Inserting into the Schrodinger 
Eq. (FT), we find that i//^i satisfies the inhomogeneous equation 

D 

HWk,\ = ^iWkfi- (42) 

We can solve for \j/ k i in terms of two independent solutions to the homogeneous 
equation, which we take as the even solution f(x) = e~ u ( x >' 2 and the odd solution 
g(x), given by 

g(x)=f(x) T-^- (43) 

Jo af z {y) 
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with the Wronskian 



Wr[f,g]=fg'-f'g=l/a. 



(44) 



It is easily verified that g also solves the homogeneous equation: Hg = 0. 
We can now construct the even solution for the first-order correction: 



Wk,i{x) 



ds 





Nk 
a 



f{x)g{s)-f(s)g(x) D 
Wr[f,g] a 2 



Wkfi( s ) 



f(x) / g{ S )f(s)ds-g(x) / f z (s)ds 
o Jo 



(45) 



In principle we could add to this an arbitrary multiple of the even homogeneous 
solution, f(x), but this just amounts to changing the normalization of the zeroth order 
solution. The even solution to the finite energy problem to first order in k 2 is then 
y/k = Nk% where 



%(x)=f(x)+k 2 a 
For large x, 



f(x) / g(s)f(s)ds-g(x) / f( S )d 



f{x)^{x/a) l l 2 - a ; x>a 



(46) 



(47) 



From this follows the large x behavior of g(x), since the integrand in Eq. (43 i is 
growing with y and the integral is dominated by the large-y contribution, giving 

{x/a) a+l l 2 



,'(*)■ 



2a 



i>a. 



(48) 



4.1 1 <a<2(l <U Q <3) 

As this point, we have to distinguish between the cases a > 1 and a < 1, as we 
did with our solvable model. We first consider 3 > a > 1, with the upper bound 
ensuring that the second moment of the equilibrium solution diverges. The calculation 
following exactly along the lines of the a > 1 calculation of the previous section. 



From Eq. ( 46 1, using the large x asymptotics of / and g, Eqs. ( 47 1 and ( 48 1, we have 



for the large x behavior of f^: 



%(x) 



l/2-a 



a > 



-k 2 a 



- (ka) 7 







l/2-a 



S {x/a) a+l l 2 

o 2aa 2a 



r f 2 (s) d s-ra) 

Jx \a/ 

(x/a) 5 /2-« (x/a) 5 l 2 - a (x/a) a+l l 2 Z 



l-2a 



ds 



4a 



4(a-l) 



2a 



2 a 



(49) 



The dominant correction term for large x scales as (ka) 2 (x/a) a+1 > 2 , since a < 3. 
Comparing this to the first, leading order, term, we see that the zeroth order term 
is dominant as long as x <$^ k~ x ' a a l ~ l ' a which for a > 1 is a more restrictive con- 
dition than x <§; \/k. Nevertheless, there is still an overlap region between the two 
approximations, Eqs. (|3~8"]l and (|49]i where we can perform the matching and find the 
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coefficients A^ and B^. The first two correction terms in Eq. (49 1 just give £?{k 2 ) 
corrections to Bk- The last term however determines A*, giving 

A k = ^(ka)- a ; st = -- LW£. (50) 



Using our value of Z for our solvable model of section III, Eq, ( 14 1, we again find 
agreement with our previous results. 

The normalization Nk is the same as before, Eq. ( 25 1, and we get for x ^s> a, 

^F a .-U(x)/2 I'™ Jk(ka\ a -Vl„(kr\ P - Dkl < 



W(x,t) « W„(x) + ^ fl e - f/ W/ 2 rdk(ka) a - l J a (kx)e 
2s& Jo 

"(a,z< 

r(a) 



e -u(x)/2 , . i/ 2 _« / Yin z 2 )\ 



For large x, W e<? (x) sw (x/a) 1 2 ^/Z, and so 



zr(o) l ' ;~*M*J r(a) 



W M « V ;; 7 ^, T(a,z 2 ) « ^(x)'-^ . (52) 



This is exactly what we found using our scaling ansatz, 

W{x 1 t)^-{ADt/a 2 ) { l 2 - a ^{z) 
a 

in our previous work 11171 . which represents the large x, large t behavior with 



(53) 



* W =^-*"IW). (54) 

This leads us to introduce the concept of the Infinite Covariant Density, which we 
define as Wjco{z,t) — W(x,t)dx/dz, giving 

W lCD (z,t) « (4Dt/a 2 ) l - a ^(z) . (55) 

Wicd, as we shall presently see, allows us to compute, to leading order, the expec- 
tation values of the higher moments of z. In fact, Wicd is the probability density 
function for z, except in a small region near z = 0, where scaling breaks down at any 
finite time. Wicd plays a dual role to W eq when considering the moments (at large t). 
The low-order moments of Wicd diverge, (including the zeroth moment, rendering 
Wicd non-normalizable), since Wicd ~ z'~ 2a as z approaches 0, but the correspond- 
ing expectation values are given correctly by the moments of W eq . Alternatively, the 
high-order moments of W eq diverge due to its algebraic decay with x, while the cor- 
responding expectation values grow as a power-law in time and are given correctly, 
to leading order, by the high-order moments of Wicd, which are finite. (The small x 
contribution, which is not correctly described by Wicd, makes an G(X) contribution, 
which is negligible to leading-order). Which moments are given correctly by Wicd 
and which by W eq is determined by Uq (i.e., a). The expectation value (\x\ q ) is given 
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by the qth moment of W eq (x) for Uq > 1 + q, (a > 1 + q/2) and by the qth moment 
of W,cd for U < 1 + q, (a < 1 + q/2): 

(2£dxx«W eq (x) U >l+q 

^ * { (4Dt)^2f"dzzmM = ««+' (f ) ^ (1+ Xff (a)z l/o < 1 +8 

(56) 
We will presently obtain the correction terms to this equation for the case q = 2. It is 



important to note that Eq. ( 56 1 implies that our system for this range of parameters 



exhibits the widely studied phenomenon of anomalous diffusion I4l l25l . 
For x's of order 1, t 3> 1, so that z <C 1, we have from Eq. ( |5"T) l, 



W(x,t)&W eq {x) + 
= W eq (x) 1 



,-U(x)/2 

2s^ 2 JQ 
a(4Dt/a 2 ) l - a n 

+ Z(a-l)r(a) 



dk{kaY a - i f{x) 



-Dk l t 



(57) 



Thus, in the central region, the solution is given to leading order by the equilibrium 
distribution, with a correction which decays algebraically in time. This implies in turn 
that a uniform approximation to W is simply 



W unif {x,t) ttW eq {x) 



r(a) 



(58) 



since the ratio of F's approaches unity as z — > 0. The F function cuts off the equi- 
librium distribution as z becomes of order unity. The physics is that the long tail of 
the equilibrium distribution takes time to establish itself, so that beyond the diffusion 
scale of \/Dt, it has a Gaussian cutoff. (For similar cases where tails take time to 
establish themselves, and so are cut off at any finite time at some large x, see lfT8l 

EH.) 

4.1.1 First-Order Correction Terms 



Turning now to the correction term in Eq. ( 57 I, it is interesting to note that the form 
of this term could be anticipated, due to the fact that the uniform solution displayed 
above is not normalized; the presence of the cutoff removes some of the probability 



from the tail of the otherwise normalized W eq : 



M r°° i / x \ i-2a r(a z 2 ) 

W{x,t)dx rj2/ W eq (x)dx + 2 / - ( - ] 1/ ; ' dx 



JO 

1- 

1 



A Z \a 
2-2a 



r(a) 



a I A \ a 

Z(a-l) \a) + zr(a)(l-a) 

2-2a 



r(a,z 2 ) 



a (A 

Z(a-l) V"^" 

a f ADt \ 

zr(a)(a-l) [~a 2 ~ ) 



a 



zr(a)(l-a) 



-r(a) 



A 



2-2a 



2-2a 



JC\2-2a x 



7 2(a-l)-z 



4Dt 



A \fl 

l-a 



2Dt 



e z dx 



\-a 



(59) 
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In the above, we broke up the integral into two pieces at x = A , with a <C A <C \f\Dt. 
If our goal is restore normalization, we can accoomplish this by rescaling W by an 
appropriate time-dependent factor which approaches unity as t — » °°: 



W(X,t)~ z ^W eq {X) r{a) , 



Z(t)=Z/ 



2\l-a 



1 



a(Wt/a 2 ) 
Z{a-l)r{a) 



(60) 



and in particular, W(0,t) = l/Z(t). 

In effect, we have added a correction term to W(x,t), This new term satisfies the 
Fokker-Planck equation in the central region to leading order in 1 /t, since the only 
violation is of order \/t relative to the correction. This is agrees of course with our 
explicit calculation in the central region, where we found precisely this additional 
term, Eq. (57i. However, in the scaling region, z ~ ^(1), the correction term does 



not satisfy the Fokker-Planck equation at all. Rather, we have to solve for the correct 
scaling form of the correction, and match it to the correction in the central region. As 
this will be important for us when discussing (x 2 ), we go through the exercise. The 
key is to write W = W m if+W\, where W„„,/ is as given in Eq. (58 1, and assume that 
the correction W\ satisfies 



Wi{x,t) ~ {ADt/ayi 2 - 2 '*^)- 



*i(z) 



>z l ~ 2a for z<l, (61) 



where the choice of the time exponent and the small z behavior is necessary for the 
matching to Eq. pl\ [j Plugging this into the Fokker-Planck equation gives 



&H + 



(2a- 1 „ \ _.. (2a- 1 



6-8a ]#•] =0 



V z J \ z L 

The solution which matches onto the central region correction is 



&!( z ) = Bze- Z U{2 -a,\ + a,z 2 )\ 



B 



ar{2 



a 



{a-l)r 2 (a)Z 2 ' 



(62) 



(63) 



where U(a,b,x) is the Kummer confluent hypergeometric function [2|. To perform 
the matching, we used the small x limit of U (a,b, x), namely 



U(a,b,x) 



r(b-i) 



r(a)x b - 1 

The uniform solution correct to second order is then 

,2\ //n.\ l-« 



X< 1 



(64) 



W(x,t)^W eq (x) 



r(a,z 2 
r(a) 



BZ 



ADt 



z a+1/2 e- z2 U(2-a,l- 



a,z 



(65) 



which exhibits the correct behavior in the scaling region, as opposed to Eq. ( 60 1 



There are of course other corrections to W(x, t) arising from the deviations of the potential from a pure 
logarithm. If U(x) ~ Uo ln(x) + ff(x~ e ) as x — > °°, there is a correction of relative order t ~ e > 2 to the PDF 
in the scaling region. This is subdominant over the correction we are considering as long as e > 2a — 2. In 
particular, for the calculation of {x 2 ) which is our prime concern, the corrections to W dominate over the 
central region contribution as long as a < 3/2, as discussed later, and the correction we consider is thus 
the dominant one in this regime as long as e > 1 . 
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It is important to note that, redoing the calculation for arbitrary xo, the result is 
unchanged (the contribution of the odd eigenfunctions, while no longer zero, can be 
shown not to contribute to leading order; this contribution is calculated in the follow- 
ing section). Thus, the solution we have found is independent of initial conditions, 
as long as the initial conditions have compact support. This is actually too strict, and 
a fast enough decay of the initial condition is sufficient to ensure universality. How- 
ever, if the initial condition exhibits a slow enough decay, then there are are always 
contributions from the region x > \j4Dt, no matter how large t is, and so universality 
is broken. For a discussion of some of the novel effects this can induce, see Ref. [ 16 1. 

In the top panel of Fig. [T] we show W(x, t) for various times for the case of the 
potential U{x) — ln(l + x 2 ), so that Uq = 2, a = 1, OC = 3/2. Together with this we 
show the equilibrium distribution, marked by t = °°. The cutoff of the distribution 
at a point which grows with time is apparent. The scaling collapse, where we plot 
t a ~ l > 2 W(x,t) as a function of x/\/4Dt is shown in the bottom panel, along with the 



infinite covariant density, Wicd, as given in Eq. (55 i, denoted by t = °°. As time 
increases, the solution in the scaled coordinate approaches Wicd- For any finite long 
time t , expected deviations (which are characterized via our uniform approximation 
Eq. ( |58] l) from the infinite covariant solution are found for small values of z- These 
deviations become small at t ^> 1; however they are important since they indicate 



that the pathological divergence of the scaling form, Eq. (55 1, at the origin is slowly 
approached but never actually reached; namely, the solution is of course normalizable 
for finite measurement times. In Fig. [2] we show W(0,t) together with the theoretical 
long-time prediction Z(t ). 

As noted above, the leading-order result for (x 2 ) is determined by Wicd- When 
comparing to numerics, one needs to take into account the highest-order correction 
term, unless the measurement time t is extremely long. There are two main correc- 
tions to the leading-order result. The first is a result of the correction term discussed 
above, which leads to a relative correction of order t l ~ a . The second is the result of 
the finite contributions for the central region, which have a relative weight of t a ~ 2 . 
Thus, for a < 3/2, the first is dominant, so that 



2\ _ / Wt \ " , r, / •*"» \ / ,2 



4Dt\ 2 - a a 3 2 f^L\ 3/2 



(x z ) « ( — ) — — „^^ w +2( —y- ) J dxx L Bze~ z U(2-a,l + a, 



a- 



) (2-a)r(a)Z \ a 2 J 



(4Dty~ a a 3 (4Dt\ 3 - 2a Df - 



2- a 


fl 3 


(2- 

2- a 


-a)r{a)Z 
a 3 


(2- 
2-a 


-a)r(a)Z 
a 3 


(2- 

2-a 


-a)r(a)Z 

a 3 



-la 



_2 



\a 2 } 



Ba 3 / dsse~ s U(2-a, l + a,s) 
o 



4Dt 

"a 27 J (2-a)r(a)Z ' V a 2 J "" r(4-2a) 

4Dt\ 2 - a a 3 (4Dt\ 3 - 2a a 4 r 2 (2-a) 

a 2 ) (2~a)r(a)Z + { a 2 ) (a-l)r 2 (a)r(4-2a)Z 2 ' ( } 

where we have used Eq. 7.621.6 of Ref. flBI to integrate the Kummer function. 

On the other hand, for 3 /2 < a < 2 the leading correction comes from the differ- 
ence between W eq and the scaling solution in the central region: 

(.wf4PA 2 -« i _ 2 r dxx 2(L(*_y- 2a _ w{x) ) . (67) 

v ' \a 2 J (2-a)r(a)Z Jo \Z\aJ qK ' ) 
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For example, for the case U(x) = t/o/21n(l +x 2 ), Uq > 2, we have Z = y/nr(a 
l)/r(a-l/2),and 



(x 2 ) a a 3 



4Df 



2-a 



1 



2 2 J (2-a)r(a)z z 



4£>A 



fdxSfx^-il+x 2 ) 1 / 2 -") 



2-a 



a 2 J (2-a)r(a)Z Z 
'4DA 2 ~ a 1 

1 "a^J (2-a)r(a)Z Z 

2-a 



2a -1 
2a -1 ^ 



f°° 2 f l 

Jo Jo 



2\-l/2-a 



4r(a + l/2) 



a" 



4Df 



1 v^(a-i) 

(2 - a)r(a)Z ~ 2(2 - a)zr{a - 1/2) 



/4DA 



2-a 



1 



1 



V a 2 J (2-a)r(a)Z 2(2- a) 
4.7.2 Numerical Tests 



(68) 



In Fig.pl we present (x 2 ) as a function of time, for the cases U(x) = f/"o/21n(l +x 2 ) 
with f/ = 5/2, (a = 7/4), and U = 3/2, (a = 5/4). We also show the first-order 
corrected predictions, Eqs. ( |68| and ( |66| , respectively, along with the zeroth-order re- 
sults. We see in both cases that while the leading-order result is quite inaccurate, due 
to the slow convergence, the first-order corrected predictions yield excellent agree- 
ment. 



4.2 a < 1 (Uq < 1) 

The calculation for a < 1 also follows along the same lines as the special case treated 



in Sec. III. Returning to Eq. (46 1, we now have 



*w-GD 



l/2-O 2 

+ k l a 



l/2-a 



o 2aa 



-ds- 



(x/fl) a+1 /2 / , 



2a 



f(s)ds 



(69) 



The last integral is dominated by its large x behavior, and so grows like x 2 ~ 2a , with a 
finite, nonuniversal contribution from the finite region. This finite contribution deter- 
mines A/, ~ k 2 ~ a , which shows that A* is subdominant to B/, for small k for a < 1. At 
this point, the calculation proceeds exactly as in the solvable example, and in fact for 



large x the result is exactly that obtained previously for a < 1, Eq. ( 36 1. The result is 



again of scaling form, with a different power-law in time as a prefactor, namely 

1 



W(x,t) « -(4Dt/a 2 )- l ^(x/Vt) 



with 



Sf(z) = 



1 



\-2a-z l 



(70) 



(71) 
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Fig. 1 Top panel: W(x,t) for various times for the potential U(x) = ln(l +x 2 ), (Uq = 2, a = 3/2, Z = 71) 
as calculated by direct integration of the Fokker-Planck equation. The convergence of the central region 
to the equilibrium distribution, W e q(x), denoted by the line t = °° is clear, along with the Gaussian cutoff 
for x larger than some threshold that moves to larger values with t. Bottom panel: The scaling collapse 
t~ , ' 2+a W(x,t) as a function of xj\j4Dt showing the approach to the infinite covariant density, Wj^q, 
as given in Eq. (55) , denoted by the line f = °°. The breakdown of the scaling for x's of order 1, which 
constitutes an ever- shrinking portion of the graph is apparent. 



Here, as opposed to the case of Wjco(z,t), ^(z) is a true normalized PDF for the 
scaling variable z. This is connected to the fact that the expectation values of \z\ q 
are time-independent for a < 1, scaling like normal diffusion. Although §f(z) is sin- 
gular at z = 0, now the singularity is integrable. Again, this singularity points to a 
breakdown of scaling for small z, here it is much less dangerous and does not lead to 
anomalous behavior. 
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• Simulation 
— Asymptotic 




•■■■»■■■«» 



400 800 1200 1600 2000 2400 

ADt 



Fig. 2 ZW(0,t) as a function of time for the case depicted in Fig. 1, together with the long-time analytic 
prediction, Z/Z(t), as given in Eq. j60|. 



Simulation 
■ Leading Order 
- 1 st Order 



- Simulation 
■ ■ ■ Leading Order 
— First Order 



Fig. 3 {x 2 ) as a function of time, as measured by simulation of the Fokker-Planck equation for the potential 
U(x) = C/o/21n(l+x ), with t/o = 2.5, (left panel), andf/o = 1.5, (right panel). Also shown are the leading- 
order power-law result, Eq. |56l, and the respective first-order results, Eq. J68) and Eq. {66) . Note that for 
small t, (x 2 ) ss 2Dt, as in simple diffusion 1171 . 



For x of order 1 , we have 

W(x,t) 



e -v{x) / 4Dt 
r(l-a)a V~o2~ 



«-i 



(72) 



which matches onto the large-x result, Eq. ( [36} and shows that the central region 
decays away slowly in time. A uniform approximation is then 



W(x,t). 



,-U(x) 



ADt 



a-\ 



(73) 



r(\ — a)a \ a 2 J 
These results of course reproduce what was found ifTTll using the appropriate scaling 



ansatz, Eq. (70 1. 



In the top panel of Fig. ffl we present (ADt) 1 a W(x,t ) as a function of x for the 
potential U(x) = !/4ln(l +ic), so that a = 3 /4, a = 1. Also shown is the asymptotic 
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result e _£/ W/r(l - a) = (1 +x 2 )- l / 4 /T(i/4), denoted as t = °°. The collapse of the 
PDF to the nonnormalizable Boltzmann distribution in the central region is clear, 
along with the cutoff at x ~ &((Dt) l i 2 ). In the bottom panel, the scaling collapse 
(4Df) 1 ' 2 W(x,f) as a function of z is shown. This collapse of course breaks down for 
small z, i-e, for x of order unity. 

As in the case of a > 1 , there is a algebraically decaying correction to the normal- 
ization, coming from the breakdown of the scaling solution at small z, which itself is 
decaying in time. If we define Z(t) = f™ m dxe~ u ( x >e~ z , we find 



z(t) 



2 / dx(x/a) 
Jo 



l -2a, 



2 dx 
Jo 



= r(l-a)a 



4Dt\ 



l-a 



-2 dx 

Jo 



\x/a) l - 2a -e- u ^ 



(74) 



The second integral of course depends on the specific form of U(x). For example, if 

U(x) = £/o/21n(l +x 2 ), we get 



dx 



(x/a) 



l-2a 



(i+(V«r) 



2\l/2-a 



(a -1/2) 



dx / ds (s + (x/a) 
(I Jo 

l 



-l/2-a 



a(a- 1/2) I ds -^M- s- a 
y ' 'Jo 2T(a + l/2) 

a^Tir(a) 



Thus, for this family of potentials, 

Z(t)=r(\-a)a(^- 



\-a 



2T(a-l/2)(l-a) 



(75) 



(l-a)r(a-i/2) 



r(a) 



(76) 



In Fig.B] we show W(0,t) as a function of time, together with the analytic prediction 
W(0,t)nl/Z(t). 

The calculation of the second moment to leading order is straightforward ifTTl : 



r(l -a)a Jo 
4D(l-a)t . 



dxx 2 z l - 2a 



ADt 



4/2 



(77) 



so that in this regime the growth of the second moment is pure diffusive in nature, 
with a renormalized diffusion constant. 

As in the case a > 1, the corrections to this result can be relatively large, as they 
decay quite slowly in time. This is apparent from Eq. f76} , where the correction term 
is of relative order t a ~ ' . While in the central region, the correction is just a change in 
normalization, in the outer region, the first-order correction has a different functional 
form. We assume the correction to W is of the form 



W 1 (x,f)«(4Df/ fl 2 ) a - 3 / 2 ^i(z); 



, 



>z l ~ 2a for z<l (78) 
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which implies that Sfi satisfies 

'2a -1 

+ ' 



it K- 



2a -1 



-4a-6 







(79) 



The solution which matches onto the central region correction induced by the correc- 
tion to Z(t), Eq. d74J>, is 



^i{z)=Cze- z U(2a- 1,1 + a,z 2 ); 
where 



C = 2J „ T ^,: X \, , (80) 



r(a)r 2 (l-a)a 2 



J = \ dx 

h 



(x/a) 



\-2a _-U(x) 



(81) 



W\ induces a correction in (x 2 
(x 2 ) w4D(l-a> + 2 
= 4£>(l-a)f 



4Df 

^2~ 



a-3/2 



o 



t/xjc 2 Cze z U{2a - 1, 1 + a,z 2 ) 



4Dt \ f°° 

—, r ) Ca 3 ds se~ s U (2a -1,1 + a, s) 
a- I Jo 



= 4D(l-a)t- 
= 4D(l-a>- 



4Dt_\ a 3 r(2-a) 
a 2 ) r(i + o) 
4DA a 2^r(2a-i)r(2-a) 

I a- 



(82) 



a 2 y r(i + a)r(a)r 2 (i-a) ' 

In Figure|6J we present (x 2 ) for the case U(x) = f/ /21n(l + x 2 ), C/ = V 2 ' ( a = 3 / 4 )> 
together with the leading and first order predictions. Again, we see that the correction 
term is important even when the measurement time is quite large, e.g., 4Dt = 10 4 . 



5 Asymmetric Initial Conditions 

If we start with an asymmetric initial condition, (a delta function at x = xq, for ex- 
ample), the first moment is not guaranteed to vanish except at infinite time. For a 
"normal" potential, with a discrete spectrum, the first moment decays exponentially. 
In the pure diffusion case, Uq = 0, however, the system never forgets the initial asym- 
metry and the first moment is preserved for all time. It is interesting then to calculate 
the decay of the first moment. The first order of business is to compute the odd eigen- 
functions. The outer solution for the odd states, which we denote by (j)k(x), is again a 
linear combination of Bessel functions: 

\J a {kx)+D k J_ a {kx)]. (83) 

y a 

The inner, i.e., x ~ &{1), odd solution we write as 0* = M^k, with 

<P k =g(x)- k 2 a g(x) f X g(s)f(s)ds - f(x) f X g 2 (s)ds (84) 
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Fig. 4 (Color online) Top panel: (ADt) l ' 1> W(x, t) as a function of .v for various times for the potential 
U(x) = (7o/21n(l + x 2 ), Uq = 1/2, (a = 3/4), showing the collapse of the central region and the cutoff at 
large x. The t = °° curve denotes the infinite time, cutoff-free result, (1 +jr 2 )~'' 4 /F(l/4). Bottom panel: 
The scaling function (ADty> 2 W(x,t) as a function of z = xj\jADt for various times, showing the collapse 
for all but a small region at small z, which disappears in the long-time limit, giving the scaling function 
S?(z),Eq. fn}. 



in analogy with Eq. (46 1. The leading order behavior of <% for large x matches with 



the small kx behavior of J a (kx) in our matching region 1 <$ix<^. l/k,so that 



C k w (fea)""^; 



2"- 1 r(a) 



(85) 



The leading order behavior of D^ comes from the constant piece in the large x asymp- 
totics of the last integral, and is nonuniversal, but by matching the k dependence of 
DkJ-a to the k 2 coefficient of f(x), we get that Dk ~ k 2+a , and so Q dominates for 
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■ ■ Simulation 
— Asymptotic 



ADt 



Fig. 5 W(0, f) versus ADt for the case presented in Fig.|4Jtogether with the asymptotic approximation for 
1/Z(r), with Z(t) given by Eq. (76). 
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Fig. 6 (x ) as a function of t, as measured by simulation of the Fokker-Planck equation for the case 
presented in Fig. [4] Also shown are the leading-order power-law result, Eq. J56), and the first-order result, 
Eq. (82). 



all a. The normalization M k is then 

{M k f : 



2<tf 2 L 



(ka 



,2a+l 



(86) 



Thus, we get, for a 8(x — xq) (xo ^ 0) initial condition, the contribution of the odd 
modes: 



W om {x,t); 



; -U(.x)/2+U(x )/2 

2^2 



dk(ka) a+l <P k (x Q )<P k (x)e 



-Dk 2 t 



(87) 
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For large t, the relevant k ~ t~ 1 ' 2 , and so Icxq is arbitrarily small and <i>k{ x o) can be 
approximated by g(xo). Then, for large x, 

„-U(x)/2+U(x )/2 I \ frr r oo 
W*«- ^ 8{X0) f-j Q dk(ka)^J a (kx)e- Dk2 ' 

-u(x)/2+u(x )/2 g{xQ) fifDtyi-ali a _ e 



2e u ^» 2 g(x ) (ADty ll2 - a _ z i 



i . i &~ z • (88) 

ar(oc) \ a J 

A more pleasing way of writing this is in term of the unnormalized even and odd 
solutions of the Fokker-Planck equation, 

X(x) = e u{x) - Y(x) = - e -~ u{x) [ X dye uiy) , (89) 

a Jo 

satisfying X(0) = l,Y'(0) = l/a and Y(x) ~x/(2aa) for large x. Then 
jodd ^ 2Y(x ) (\Dt\-^- a , 



w -m^fwA-^) * ■ (90) 

This odd solution has the scaling structure W odd ~ t$ & dd{z), j3 =— 1/2— a. Whereas 
for the even solution, we could guess the exponent j3, since we could demand that it 
match the time-independent leading order solution, here there seems to be no way to 
guess j3 without using our eigenvalue expansion. Plugging the scaling ansatz into the 
Fokker-Planck equation gives the scaling equation for & dd- 

*« , f 2 "" 1 ,o,W J^*„ 2a - l \& 

2 



= &" dd + +2z ^ odd + 2 + 4a — ] & odd , (91) 



\ z J \ z 

2 

for which ,^ dd{z) = ze~ z is clearly a solution. For x's of order 1 

-U(x)/2+U(x )/2 



W odd (x,t) = Wl j Q g(x Q )g(x)(ka) 2a+l e- Dk ' 



-l-a 



_ 4aZY Q (xo)Y (x) (ADt\ 

X(x Q )r(a) {a 2 J ' ^ 

so that for 1 <C x <C \/4Dt, this approaches our previous result. A uniform approxi- 
mation for the W odd is then 

4aY (xo)Y (x) /WV 1 -^ 
aX(x )r(a) \a 2 J K ' 

We see that for all a this decays in time as a power law which is faster than the decay 
of the even part. Furthermore, it is xq dependent, as opposed to the even part which 
to leading order is independent of xq. 
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Fig. 7 (x) versus 4£>f for the potential U(x) = l/4ln(l +x 2 ), for the cases x — 1, (blue dotted line), 10, 
(red long dashed line), and 40, (purple short dashed line), together with the asymptotic approximations, 
Eq. (94), shown as the corresponding solid lines. 



The calculation of the first moment is then straightforward, since it is dominated 
by the large x behavior: 



4Y(x Q ) 



a- 



aX(x )r(a 

a^/7zY(x Q ) I -\:I)i 
X(x )r{a] 



<?) 



-1/2-a 



dxxze 



1/2- a 



(94) 



In particular, when Uq = 0, i.e., a = 1/2, the first moment does not decay at all. 
Since Y(xq) vanishes at xq — and grows for large xq as xq, the first moment grows 
as jcq", so starting with a fairly large xq should make measurement of this decay much 
easier. We see this in Fig. FT] where we present (x) as a function of time, for xq = 1, 
10, and 40, for the potential U(x) = i/4ln(l +x 2 ). We see that even though the time 
for the asymptotic decay to set in gets much larger as xq increases, the value of the 
first moment at the onset of the asymptotic scaling regime is still much larger for the 
larger xq. 



6 Conclusions 



We have explicitly constructed the probability density function for diffusing Brow- 
nian particles in an attractive logarithmic -type potential, V(x) ~ Vblnjc, Vq > 0. The 
system exhibits different modes of behavior as the strength of the potential is varied. 
For very weak potentials, Uq = Vo/ksT < 1, the particles diffuse as t 1 ' 2 (to leading 
order) with a re normalized diffusion constant, and with algebraic corrections. The 
Boltzmann distribution in this case is nonnormalizable, yet nevertheless the PDF, for 
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x <C \/ADt, is prorportional to the Boltzmann density, with a algebraically decaying 
pre factor. For 1 < Uq < 3, we see that both the normalizable Boltzmann distribution 
and the infinite covariant density are necessary for a full description of the leading- 
order dynamics. The probability density approaches the equilibrium distribution for 
x's smaller than 0{\fDt) in magnitude. For larger x's, the PDF decays as a Gaussian. 
The second, and higher, moments of the distribution are sensitive to the cutoff, grow 
algebraically with time, and can be calculated from the infinite covariant density. For 
larger Uq, the second moment saturates at its equilibrium value, but higher moments 
are sensitive to the cutoff and grow as a power-law in time. These cutoff-sensitive 
moments are calculable via the infinite covariant density, which is the scaling limit 
of the PDF. Thus, in all cases, the Boltzmann distribution describes the x <C \/4-Dt 
regime, and a scaling function describes the PDF at larger x. The leading-order calcu- 
lation of the moments is subject to corrections which decay quite slowly in time (i.e. 
with a power less than one-half), and so for reasonably large times, the leading-order 
prediction is quantitatively poor. Upon adding the first-order correction, the predic- 
tion of the Oth and second moments is seen to be quantitatively accurate. The first 
moment of the distribution, unlike for pure diffusion where it stays constant, decays, 
albeit very slowly, as t~ u °> 2 . 
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